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Abstract 

A multi-channel scattering problem is studied from a point of view of integral equa- 
tions system. The system appears while natural one-particle wave function equation 
of the electron under action of a potential with non-intersecting ranges is considered. 
Spherical functions basis expansion of the potentials introduces partial amplitudes and 
corresponding radial functions. The approach is generalized to multi-channel case by a 
matrix formulation in which a state vector component is associated with a scattering 
channel. The zero-range potentials naturally enter the scheme when the class of operators 
of multiplication is widen to distributions. Oscillations and rotations are incorporated 
into the scheme. 



1 Introduction 

In 1936, Fermi proposed zero range potential (ZRP) model to study neutron scattering 
in hydrogen- containing substances. Since then, ZRP approach have been developed to widen 
limits of this pioneering treatment (for a review see Demkov and Ostrovsky 0, Drukarev 
Albeverio et al. The advantage of the theory is the possibility of obtaining an exact 

solution of scattering problem. Recently, Baltenkov [|10] have generalized the ZRP method for 
the case of non-zero orbital moments (see ^ )including combinations of potentials . De Prunele 
IP have proposed other solvable non-zero range potentials, which involve higher partial waves. 

The aim of this paper relates to the other limitation of the real scattering phenomenon. 
Following the main ideas of the ZRP method we introduce a matrix ZRP potential, (Section 

for a multichannel problem (see, for instance. Lane The matrix ZRP is conventionally 



represented as the boundary condition on the matrix wavefunction at some point. Presented 
matrix potential generahzes ordinary matrix ZRP, which was proposed by Demkov and Ostro- 
vsky [|TT| (see also []T^), for the case of the s,p,d, etc. target states. We also present simple 
method of deriving scattering amplitude for a multi-center problem (Section and consider 
two matrix ZRPs problem (Section^). As an important example we consider applications to 
a diatomic molecule. The formulae for differential and integral cross sections of the electron- 
vibrational excitations are summarized in the Section |^. In the Section ^ we present the results 
of our numerical calculations for the molecule H2 as well as plots for electron-vibrational cross- 
sections. 



2 Nonoverlaping potentials and multi-center problem 
2.1 One-center case 

Let us consider the scattering problem for matrix wavefunction y{r, Kiiq) and short-range 
matrix operator, which is U ai r < d and zero at r > d. The atomic units are used throughout 
the present paper. In the interior (i.e. at r < d) 

iK\r—r'\ 

y{v, Kno) = e^^"-^ - i2n)-' [ — f/y(r', Kno)dv', (2.1) 

J |r— r 



where K = diag(/c„), kn = yk^ — 2En are electron momenta in the channels n, En are excitation 
energies, and no indicates incoming electron direction. The matrix wavefunction y{r, Kuq) in 
the exterior region (i.e. at r > rf) is a solution of the Helmholtz equation and is determined by 
smooth matching condition at the boundary (i.e. at r = d). 

In order to construct solutions in the exterior region we introduce a matrix function H{r, u; K) 
as a solution of the equation 



gi/C|r-r'| 



|r-r'| Jn 



H{v, u; K)e~'^'' ''' dn^, r' < r, (2.2) 



where u is unit vector. This diagonal matrix function is invariant under rotation transforma- 
tions of the vectors r, u and has the following asymptotic behavior at infinity 

H(r,u;K)'^°^^6(n-u), (2.3) 
r 

here n = r~^r, and 6{n) denotes delta function of the angle variable n. Solution in the exterior 
region, i.e. y{r, Kuq) at r > d, can be constructed in terms of the matrix function H{r, u; K) 

y{r, Kno) = e^^""-"" + / H{r, n; K)F{u, no; K)dn^, (2.4) 

Jn 

where matrix amplitude F{n, uq] K) is given by the expression 

F(n, no; K) = -{27i)-' J e-^^"-'-f/y(r, Kno)dr, (2.5) 
in which n is outgoing electron direction. 
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2.2 Multi-center case 



Consider nonoverlapping short-range matrix operators e"*^ "^* Ui e^^'^\ here p describes 
the momentum operator —iV. Suppose that the interaction regions are hmited by spheres of 
the radiuses di with centers at the points Rj. Requirement of nonoverlapping is written down 
as di + dj < |Rj — Rj|. Denote the matrix wavefunction in the region |r — Rj| < di by the 
expression e^*P''^'\&i(r). The matrix wavefunctions \l/i(r) satisfy the following equations in the 
region r < di 

N i_R'|r-r'+R,;-Rj| 

^.(r) = e^^""-(^+^') - (27r)-i ^ / U,^j{r') dr', 

j=l |r— r'+R,;— Rj 

where integration must be performed over r' < dj. Taking into account the condition of 
nonoverlapping and the equation ( p.2|) we obtain the expressions 

iXIr-r'l 
|r-r'| 



iK\r—r'\ 

^.(r) = e^^"« ('-+^') - (27r)-i / f/,^i(r')rfr'- 

J |r— r' 



N 



(277)-^^ / / H{Ri - Rj,u; K)e'^''<'-'''^Uj^j{r')dr'dn^. 

To express the matrix wavefunctions ^'i(r), let us use solutions of the correspondent one-center 
problems. In terms of the matrices yi{r, Ku), the matrices ^'i(r) become 

^,(r)= / / y,{r,Kn)F-\n,n';K)Q{n',no)dnndnn', 
Jn Jn 

where Cj(n',no) are some constant matrices. Taking into account Eqs. (|2.4| ) and ( p. 51) we 
obtain 

N 

/ F-l(n,n';ir)a(n^no)rffi„. = e^^"''-^'(5(n-no) + E^(R.-R.,n;ir)C,(n,no). (2.6) 

In order to calculate multi-center F-matrix for the nonoverlapping potentials we only need to 
solve this matrix equation in N unknown matrices Cj(n, Hq). The multi-center F-matrix is 
given by 

N 

F(n, no; i^) = E e-^^"-^'a(n, no). (2.7) 

1=1 

We will omit from now the argument K of the matrix F-function. 

If we take into account only s-wave point interaction then Eqs. (pl6|), ( p.7| ) are reduced 
to usual equations of the ZRP theory. Another solvable model was recently proposed by de 
Prunele |Q . In this model nonoverlapping separable interactions are localized on the spheres of 
the radiuses di with centers at the points Rj. It is clear that the most essential aspects of this 
model are obtained also by using the Eq. ( |2.6| ). Thus, presented formalism can be a basis of a 
new solvable models creation. 
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3 Matrix ZRP and boundary condition 



Now we examine closely the electron-center interaction on the multi-channel level of approxi- 
mation. If interaction preserves a total angular momentum of electron-center system then the 
matrix partial waves are also eigenfunctions of the total and projection angular momentum 
operators. Further we consider for simplicity a case of the zero total angular momentum. As- 
sume Ini^n denote angular momentum quantum numbers of isolated center in the channel n. 
The isolated center states may be characterized as d, etc. for 1^ = 0, 1, 2, etc. respectively. 
Hence matrix function y{r, Kuq) ~ ^(n), where A{n) is the diagonal matrix 

A(n) = diag(F,„„„(n)). (3.8) 

Thus, /„ is also orbital momentum of the incident electron in the channel n; rrin is its projection 
on the axis z. 

Following the main idea of the ZRP method, we suppose that matrix interaction U is 
localized at the center point (the point can be considered also as a sphere of radius d 
0+). Therefore the matrix y{r,Kno) is a superposition of the regular and irregular Helmholtz 
solutions. The asymptotic behavior at zero is given by (further we use for simplicity the 
notations like diag(r^") = r^^^^i''") and diag(?/^") = diag(|/n)*^^^S^^"-'.) 

y{r, Kno) ~° A{n) (^{2L - E)\\ r"^"^ - J^J^^^W^ A+(no), (3.9) 

where E is unit matrix, L = diag(/„). The Hermitian matrix W fixes relation between any 
regular and irregular Helmholtz solutions. The reactance matrix can be expressed in the terms 
of the matrix W as —K^'^^^'^W'^K^'^^^'^. In the Eq. ( ^lOD we hold only the leading terms of 
regular and irregular solutions. Such asymptotic leads to the boundary condition on matrix 
wavef unction 



\ 2L+E 



dr I 
~2^L\W 



A^{n)y{v,Kno)dnr, 



(3.10) 



r 



r=0 

L+E 



{2L + E)\ 



A+{n)y{r,Kno)dn^ 



r=0 



Breit proved (p[, see also H, 0) that the ZRP and hence matrix ZRP can be also introduced 
as pseudopotential. 

Imposing the boundary condition on the integral representation of matrix function y{r, Kuq) 
(see Eq. (|2.4| )) we derive the matrix amplitude 

F(n, no) = 47r A(n)FA+(no), (3.11) 

where E is given by 



'^WqK'^ + iK) \ (3.12) 



and Wo = {—i)^Wi^ is also hermitian matrix. 

Let us consider special cases. 
• One-channel ZRP 

For the case of any state {Iq = I) and Wq = a scattering amplitude becomes 



a + ik'^^'^^ 
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where a = 1/a - inverse scattering length. The expression coincides with scattering amphtude 



F, which can be calculated for scattering by GZRP ||T0[. Assuming Z = we get scattering 
amphtude of isolated ZRP ||^. 
• Two-channel matrix ZRP 

Let us consider ground s {Iq = 0) and some (/i = /) excited states. Equation (|3.12|) goes into 
the inverse matrix amplitude 

L=(^, ? V F-^ = -f «o + fo V (3.13) 

where k^^ki are related by the energy conservation low and q;o,q;i,c are real parameters. In 
the case / = we obtain one-center matrix amplitude for two s states 0. 

4 Two matrix ZRP problem 

In order to study thoroughly the multiple scattering on a multi-channel level of approximation 
we consider two matrix ZRP problem. Assume two-center matrix potential U is localized on 
the points Ri = R, R2 = — R and satisfies the following parity requirement 

Pf/p-i = Sf/S"\ (4.14) 

here P denotes inverse operator Pf{r) = /(— r), and S = diag(?7„) is the matrix of parities. 
Let us represent U as 

where Ui, U2 are matrix ZRPs, and 2R is a distance between the centers. This representation 
leads to 

U2 = T.PUiP-^T.. 
Taking into consideration these expressions, we conclude that 

Ffi(n,no) = F-i(n,no), (4.15) 

F2-i(n,no) = SF-i(-n,-no)S. (4.16) 

In the body-frame, where polar axis z is taken along R, matrix -F(n, no) can be represented by 
Eq. (|3TlD . Denoting 

Ci(n,no) = (47r)i/2^(n)Ci(no), ^2(11, no) = (47r)i/2 A(-n)EC2(no) (4.17) 
we obtain for matrixes Ci(no), C2(no) the equations 

r F-iCi(no) - SifC2(no) = (47r) VM+(no)e^^"°-^, 

\ F-iC2(no) - SffCi(no) = [A-nf^ SA+(-no)e-^^"o-^, ^ ' ' 

where matrix H (argument KR is omitted) is given by 

H = An [ A+(n)H(2R, n; K)A(-n)dnn. (4.19) 
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In order to calculate matrix H we can use the following expansion 

K °° 

H{2R, n;K) = — Y] ^^(2A + l)hx{2KR) Px{R~^R ■ n), (4.20) 
4^ A=0 

/ 1 (i \ e"' 

where h\(x) = { — — are Riccati-Hankel functions, P\{x) are Legendre polynomials. 

y X ax J X 

The integral ( 4.11^ ) can be evaluated in terms of 3-j symbols or Clebsch-Gordan coefficients. 
Thus, determination of the matrixes Ci(no), C2(no) by using Eqs. ( [4.18|) reduces to solving a 
linear system of algebraic equations. The final result for the matrix amplitude -F(n, no) becomes 

(27r)-iF(n,no) = 

(A(n)e-^-^"-^ + EA(-n)e^^"-^)(F-i - S/7)-i(A+(no)e^^"«-^ + SA+(-no)e-^-^"«-^) + 

(A(n)e-^^"-^ - SA(-n)e^^"-^)(F-i + S/7)-i(A+(no)e^^"o-^ - SA+(-no)e-^^"«-^). 

(4.21) 

For purposes of illustration, we represent some examples. 
• One-state level of approximation 

Two-center electron-molecular interaction may be approximated by two ZRPs. The simplest 
choice for E-state is two s-centers, i.e. ^(n) = (Air)^^^'^. If one-center inverse scattering 
amplitude and H are 

F-^ = -a-ik, H = , 

then two-center scattering amplitude is given by the expression 

cos(A;n-R) cos(A;no-R) sin(A;n-R) sin(A;no-R) 
F(n, no) = -2 ^rj^ 2 ^rj^ . 

a + ik -\ — — — a + ik 



2R 2R 

The result coincides completely with scattering amplitude for two ZRPs 
• Two-state level of approximation 

The problem of electron-impact excitation of the S, 11, A, etc. molecular states can be con- 
sidered in two-state level of approximation. Suppose, for a simplicity sake, molecule have even 
ground S-state (i.e. tjq = 1) and any excited state, which may be S, 11, A, etc. (i.e. m = 0, 1, 2, 
etc. and r^i = 1 or —1). Let us approximate two-center molecular interaction by two matrix 
ZRPs. One-center matrix F is given by Eq. (|3.13|) , where / > m. We introduce the following 
notation for short 



9o{x) = ao + iko + X 



^^21 (4.22) 



e^[x) = a^ + ikf^^ + x t^'+\2X + l)kf+'hx{2k,R) j |r^„(n)pPA(n ■ e,)rffi, 

A=0,2, ''^ 

where is unit vector orientated along polar axis z. Thus, formula ( |4.21|) results in the elastic 
scattering amplitude (i.e. for a transition — 0) 

cos(fcon-R) cos(A;ono-R) sin(fcon-R) sin(fcono-R) 
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Assuming r/i = —(—1)', we obtain electron-impact excitation amplitude (for a transition 0^1) 

^ , , ^ ^7i^r I X /cosf/cin-R) cos(A;ono-R) sin(A;in-R) sin(A;ono-R) \ 

Fio n, no = 4v^cA; n \) , J , ' ^ + J , \ ^ . 4.24 

V ^o(^?o)^i(^i) - c2 ^o(-^o)^i(-^i) -c^ J 

In the case rji = (—1)' one becomes 



^ / N ^t/t^ / ^ /sinf/cin-R) cosf/coHo-R) cosf/cin-R) sinf/coHo-R) \ ..^^s 



According to general theory |P, ^, a resonances occur in the vicinity of the poles. The poles 
of the matrix amplitudes correspond to the solutions of the equations 6'o(±r7o)^i(='='7i) = c^- 
Considered as a function of the spacing on centers 2R these energies represent the adiabatic 
potential curves of the negative ions or quasistationary states. In the matrix ZRP model the 
poles may reproduce both shape and Feshbach resonances. 

The matrix amplitudes are represented in body-fixed frame. The amplitudes in other frame 
can be obtained by frame rotation. The frame-transformation is reduced to the simple sub- 
stitution Yim{n) — Y.m'=-i ^mm'^im'{T^) siuce the scalar products Hq-R, n-R are invariants of 
a frame-rotational transformation. Here -D^m' denote the rotation matrix elements (so-called 
D-functions, 0). 



5 The adiabatic approximation 



The adiabatic approximation can be used in order that to incorporate the motion of the nuclei 
into the theory. Initially this approximation was applied by Drozdov |13|, Chase [n|, and 
Oksyuk The adiabatic approximation within the framework of the ZRP model was devel- 
oped by Demkov and Ostrovsky and Drukarev and Yurova |Q. This approximation allows 
to express the electron-vibrational transition differential cross section (DCS) via the electron 
transition amplitude on the space-fixed matrix ZRPs: 



da 

dn' 



Ot;o) = (47r)-^M„^ / |(nt; |F„o(n, no, R)| 0^;o)N^^ 



(5.26) 



where n,v,0,vo represent electron and vibrational quantum numbers for final and initial states, 
{nv\ ■ \Ovo) denotes integral of the vibrational harmonics, and M„ is the orbital angular mo- 
mentum projection degeneracy factor of the final target state n. 

The integral cross section (ICS) for electron-vibrational transition is obtained by the inte- 
gration over the scattering angle: 



a(nv ^ Ovq) = [ -^{nv ^ Ovq) dQn 
Jn dil 



(5.27) 



The summation over the vibrational states (generally - including continuous spectrum) gives 
the electron transition ICS 



a[n 



Ovo) = {A7T)-'Mn^(ovo I I |F„o(n,no,R)|'t/^]A„ 
fco \ J^ 



(5.2J 



In these equations the integrals over angles can be reduced to the spherical harmonics sums, 
which are most suitable for numerical calculation. 
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• The electron-vibration transition DCS for two-channel problem 

It is convenient to represent the amplitudes Eqs. ( [4.24| , [4725D in coordinate frame, where the 
polar axis is taken along vector n. As it proved above, the rotational transformation is reduced 
to the substitution Yim(n) — * Yim(R)- Expanding the amplitudes in the series over spherical 
harmonics we obtain the expression 



da 



7, oo oo min(A,A') 

(47r)2 M^-^ J2Y: E Q>^,Qy, / \YumA^)\'Y*^in)Yy,iR) dn^, 



A=0 A'=0^=-min{A,A') 



where integral of the spherical harmonics can be evaluated in terms of the Clebsch-Gordan 
coefficients, and are given by 

vam — ^ X 

Yx^iKn + fcono) {nv \ - Z^)) jxi\Kn + kono\R)\Ovo) + 



where 



(-0i(r/i), ck[' 



z(-) 



• The electron transition ICS for two-channel problem 

Substitution Eqs. ( ^ , [41251) into Eq. ( ^1281) yields the following ICSs 



ain Of I 



(-l)'"r7„5,„„„(2A;„i?) + 



Blmix) 



(l - ^0^^^^^) (1 - (-l)Sn5.„™„(2A:„i?) 
i\'2X + l)jx{x) I m„(e)|2p,(e-e,)dne- 



6 Applications and discussion 

The integral cross section for electron-vibrational excitation of a state of H2 is plotted in 
fig. |I] for a number of values of Q;o,ai,c, i? (see [0), which are regarded as constant in the 
range of interest. 



7 Conclusion 

Among the most important aspects of the paper are the development of the matrix ZRP 
theory, the calculation of the differential and integral cross sections for the electron-vibrational 
transitions, and investigation of the matrix ZRP possibilities. The nonoverlapping condition 
essentially simplifies a scattering problem solving. Thus, the determination of the multi-center 
matrix amplitude reduces to solving of the system (see Eqs. (p.6|), (|2.7|)) of integral equations. 
The number of scattering centers defines the number of equations in the system. 
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Figure 1: Total cross sections for excitation of f = 0, 1, .., 7 vibra- 
tional levels of the a^^g state of H2. 
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